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Abstract. In young star clusters, the density can be high enough and the 
velocity dispersion low enough for stars to collide and merge with a significant 
probability. This has been suggested as a possible way to build up the high- 
mass portion of the stellar mass function and as a mechanism leading to the 
formation of one or two very massive stars (Af* > 150 M©) through a coUisional 
runaway. I quickly review the standard theory of stellar collisions, covering 
both the stellar dynamics of dense clusters and the hydrodynamics of encounters 
between stars. The conditions for collisions to take place at a significant rate are 
relatively well understood for idealised spherical cluster models without initial 
mass segregation, devoid of gas and composed of main-sequence (MS) stars. In 
this simplified situation, 2-body relaxation drives core collapse through mass 
segregation and a coUisional phase ensues if the core collapse time is shorter 
than the MS lifetime of the most massive stars initially present. The outcome of 
this phase is still highly uncertain. A more realistic situation is that of a cluster 
still containing large amounts of interstellar gas from which stars are accreting. 
As stellar masses increase, the central regions of the cluster contracts. This 
little-explored mechanism can potentially lead to very high stellar densities but 
it is likely that, except for very rich systems, the contraction is halted by few- 
body interactions before collisions set in. A complete picture, combining both 
scenarios, will need to address many uncertainties, including the role of cluster 
sub-structure, the dynamical effect of interstellar gas, non-MS stars and the 
structure and evolution of merged stars. 



1. Collision and Interaction Timescales 

We first consider the timescale for direct collisions between stars. We assume 
two populations of stars, with stellar masses mi and m2, stellar radii ri and r2, 
and Maxwellian velocity distributions of (ID) dispersions ai and a2- The mean 
time for a star of type 1 to collide with a star of type 2 is 



W = { V87rn2cri.ei(ri + r2) 



^ ^ G{mi + 7712) 



(1) 



with cr^j,2 = '^1 + ^2 ( Binnev &: Tremaind 19871 ). The second term in the square 



brackets originates in the mutual gravitational attraction of the stars. This 
gravitational focusing dominates over the purely geometrical cross section (rep- 
resented by the term 1 in brackets) unless the velocity dispersion is similar or 
higher than the escape velocity from the surface of a star, i.e. 500 — 1000 km s~^ . 
Such relative velocities are presumably only reached in the vicinity of massive 
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black holes (e.g., Freitag &: Bend 20021 . 20051 ). For other situations, the collision 
time can be written 



^cou ^ 5 Gyr 



lO^pc" 
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10 km s ^ ri + r2 rrii + m2 



(2) 



In Fig. [H I show (with dashed lines) the stellar densities required for a massive 
star to experience on average one collision during the MS. Typical values are in 
excess of lO^pc"^ or 10® pc^^ for a star of 10 Mq or 120 M©, respectively. 

Stars passing each other within a few stellar radii at low velocities can 
dissipate enoug h orbital energy in tides to form a bound binarv (e.g.. Fabian. 



Pringle, &: Rees 19751 : iKim &: Leelll999l ). Such a "tidal binarv". being born very 
tight, is likely to merg e quickly, or, at least, to enter a common-envelope phase 
(iDale Dav ied l2006h. The merger cari also be induced by the perturbation 
of a passing star ( Bonnell h Batd l2005l ). In any case, because tidal dissipa- 
tion decreases steeply with increasing closest-approach distance, formation and 
merging of tidal binary can reduce the effective value of the collision time by a 
factor ~ 3 — 5 at most. 

A more promising way to get interesting collision rates is to consider binary 
interactions. The timescale over which a binary with separation a and mass 
mbin suffers from a close interaction with a field star is simply 
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where m is the individual mass of the field stars and their number density. 
For a sufficiently hard binary, i.e. when the total orbital energy of the binary and 
impactor system is negative, most close interactions are "resonant" . This means 
that the three stars orbit each other in a chaotic fashion for several dynamical 
times within a volume of order until this meta-stable system decays into a 
binary and single star if stars were point masse (e.g.. iHeggie fc HutI [2003l 'l . 
There is in fact a high probability tha t two of the stars coll ide and merge during 
a resonant interaction. For instance, iFregeau et al.l (|2004l ) showed that, in the 
hard regime, the cross section for a collision to occur during an encounter of a 
circular, a = 1 AU binary consisting of two 1 Mq MS stars with a similar star 
is about |7raGMoz;~p corresponding to a hard sphere of radius a/4. If massive 
binaries are born with eccentricities « 0.9, as a result of accretion of interstellar 
gas, more distant, non-resonant , interactions with pas sing stars can suffice to 



cause a merger, as suggested bv lBonnell fc Batel ^M)- 

At relative velocities below about 5kms~^, encounters involving two mas- 
sive stars, one or both of them surrounded by a massive disc, can result into a 
bound binary if the closest approach distance d is smaller than a few tenth of the 
disc ra dius, which corresponds typically to d ~ 10 — 100 AU (jMoeckel &: Bally 
2007bl ). With d substituted for a, Eq. [3] suggests that this mechanism might be 
responsible for the high proportion of massive stars with a companion of similar 
mass in clusters such as the ONC, despite the sh ort life time of massive discs 
(iMoeckel Ballvl l2007al : iPfalzner QlczakI l2007l ^ . 



^Binary-binary interactions occur in a similar fashion. 
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The Outcome of Stellar Collisions 



Collisions between single MS stars have been the object of extensi v e num erical 
studies, starting with the 2-D grid simulations of ISeidl Sz CameronI (|l972( ). with 
most works based on t he Smoothed Particle Hydrodynamics (SPH) method (e.g., 
Freitag fc B"en3l2005l : see the contributions in Shara 2002 and the "MODEST" 
web page^ for more references on stellar collisions). At low velocities, typical 
of globular clusters or young clusters, any collision between MS stars lead to a 
merger with very little mass loss (< 10%, typically a few %). Neglecting the 
possibility of the formation of tidal binaries, stars can therefore be treated as 
sticky spheres to good approximation. 

In this velocity regime, one can determine how the mass elements of the 
parent stars sort themselves in the collision product, using a algorithm based on 
Archimedes' principle, to determine the internal structure of the merger prod- 
uct with go od accuracy, avo i ding t he computational cost of an hydro dynamical 



simulation ( Lombardi et al.l l2002l : iGaburov et al.ll2007bl ). In case of collisions 



between mass of unequal masses, most of the material from the lower object 
sinks to the centre of the merged star. 

One niaior remainirig uncertaintv is the evolution of collision products ( Sills 



et al. Il997l . l200ll . l2002l : iGlebbeek Pold 120071 ). When the merger has settled 
down to a new hydrodynamical equilibrium, it is significantly swollen compared 
to the size of a MS star of the same mass, and typically spinning at a rate close 
to breakup. Some efficient mechanism, such as winds or magnetic locking to a 
disc, has to operate to shed most of this angular momentum but comparatively 
little mass as the star contracts back to t he MS, lest it be nearly corn pletely 
ground down by "centrifugal evaporation" ( Sills. Adams. Sz Davies [20051 ) . 

From consideration on cross section and the fact that most massive MS stars 
are members of binaries (or higher-order systems), it seems that most collisions 
in a young cluster occur during binary interactions. Unfortunately, because 
of the added numerical complexity and the enormous parameter space to be 
explored, only very f ew researchers have carried out hydrodynamical sinr ulations 
of such interactions ( Goodman Sz Hernquist|[l99ll : iLombardi et al.ll2003l ). These 
works made it clear that multiple mergers are relatively likely during a binary 



interactions due to the increased size of collision products (see also Fregeau et al, 
20041 1. 



In a very young cluster, one has to consider collisions involving stars still 
on the pre-MS. Low-mass stars (m < 5 — 10 Mq) take much more time than 
massive ones to con tract to the MS and can therefore be significantly larger that 
massive stars (e.g.. [Pallal |2002| ) . Only few SPH simulations involving pre-MS 
stars have been published so far (|Zinnecker k Batdl2002l : iLavcock fc Sinil2005l : 
Davies et al.l l2006). They show how, in a close encounter with a massive star, a 
low-mass pre-MS star is tidally disrupted and end up forming a disc around the 
massive star. As mentioned above, such massive discs offer a very large cross 
section for the formati on of tight binaries which may later merge, a sequence of 
events investigated bv lDavies et al.l (j2006l ). 



^ "Modelling Dense Stellar Systems" ; see jhttp : / /www . manybody . org/modest/[ Follow the link 
"WG4" for the pages on collisions. 
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Figure 1. Stellar densities required for high collision rate. The dashed lines 
show the density at which the collision time icoU for a star of mass Wmax = 
10 M0 (left panel) or 120 Mq (right panel) is equal to its MS lifetime tms, at 
the centre of a non-evolving King-model cluster with Wq — 3, 8, 12 (from top 
to bottom). The dependence on cluster mass (and on Wq) comes from the 
dependence of tcoii on the velocity dispersion. The solid lines show the central 
density corresponding to a core-collapse time ice equal to tms- To estimate 
tec, we assume a Salpeter mass function from 0.2 Mq to mmax- The dots on 
the solid lines (right panel) are a rough estimate of the number of stars above 
which our estimate of tec is robust. 



3. First Scenario: Relaxation-Driven Core Collapse 



Several authors have studied how 2-body relaxation can drive the evolution of 
a cluster to a stage of very high central density, leading to stellar collisions 
( Portegies Zwart et al. 19991: Portegies Zwart &: McMillanI 20021: Portegies Zwart 



et al. I2OO4I: iGiirkan. Freitag Rasidl2004l: iFreitag. Rasio. fc Baumgardtl l2006bl: 



1 I ^ ^ I f M " ' ' '^^ I f I — — " 

Freitag. Giirkan. &: Rasid l2006a: Giirkan, Fregeau. &: Ra sio 2006: Gaburov, Gua- 
landris, &: Portegies Zwart ,2007a). Two-body relaxation is the process through 
which stars in a populous cluster (A^ ^ 10) exchange energy and angular momen- 
tum over a time scale significantly longer than their orbital time, by deflecting 
each other's trajector}!^. Roughly speaking, the relaxation time is the timescale 
for orbital parameters of a an "average star" to change completely as a result of 
a very large number of uncorrelated small-angle h yperbolic 2-body encounters 
with all other stars (e.g., Binnev Tremaine 19871 ). Its local value is given by 
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■^In a smaller group, 2-body relaxation cannot be defined in a useful way because the energy 
and angular momentum of a star are generally not conserved on a dynamical time. 
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Figure 2. Core collapse times for star clusters with a variety of mass func- 
tions. Shown are previously unpublished results of Monte-Carlo simulations 
by the author (circles) and of A^-body simulations by H. Baumgardt (crosses). 
In most cases, a Plummer model with a MF dN/dm cx m^'^-^^ was used. 



where A ~ ln(0.02A^), N is the number of stars in the cluster, a the ID velocity 
dispersion and (m) the mean stellar mass. 

In a cluster with a mass function (MF), the massive stars lose kinetic energy 
in favour of lower-mass objects through 2-body relaxation. This would eventu- 
ally lead to energy equipartition if it was not for the cluster's self-gravity. As 
they lose energy, the most massive stars concentrate to the centre, a process 
known as dynamical mass segregation and which reduces to dynamical friction 
in the limit of very large mass ratio. For any reasonable MF the central subsys- 
tem of massive stars becomes self-gravitating before equipartition is established. 
This makes thermal equilibrium impossible to achieve because the massive sub- 
system gets denser and hotter at an accelerated rat e as it g i ves u p energy to 
the rest of the cluster, an instability first predicted by Spitzeil ( 1969|) . Using the 



Mon te-Carlo stellar dynamical method invented bv Henon QOTla^, Giirkan 
et al. (|2004l ) showed that in rich clusters this mass-segregation-induced core col- 



lapse occurs on a timescale tec dictated by the initial central relaxation time 
trc(O), largely independently of the cluster structure, 

t,,^2t,M- for /.= ^>5. (5) 

This relation is shown in Fig. [2] for a set of cluster simulations. 

The solid lines in Fig. [T] show what is the minimum initial central density of 
a cluster of given mass and concentration for the core collapse time to be shorter 
than the MS lifetime tms of the most massive stars. In Monte-Carlo simulations 
where this is the case, the core collapse proceeds until the massive stars in the 
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central region start colliding, either as single stars if there are no primordial 
binaries (a rather unrealistic idealisation) or during binary interactions. When 
tec > tms, the gas lost by massive star in the post-MS evolution causes the core 
of the cluster to re-expand before any significant number of collisions can occur. 

Simulations carried out with the more accurate but very time-consuming 
direct A^-body algorithm in the regime > 40 and N < few x 10^ indicate 
a more con straiiiing condition for a collisional phase to occur fe.g.. Portegies 
Zwart et al. \2004 ]. This is probably because the type of core collapse described 



here can only occur if there are initially a sufficient number of stars with masses 
close to ?7ijnax in a region not much larger than the core. In Fig. [1] we indicate 
with dots on the solid lines the cluster masses below which there are on average 
fewer than 5 stars with a mass between O.Smmax and mmax in the core initially. 
This is not ci significcint limitcttion for iriYnax — 

10 M0. 

The outcome of the collisional phase is still a subject of debate. Dynamical 
simulations which only take into account collisional mass loss and assume that 
the merger product returns to the MS immediately after a collision show a 
runaway growth of one or two stars to masses nij-g, of a few hundreds or thousands 
Mq. Such "very massive stars" (VMSs) have been suggested as progenitors 
of intermediate-mass black holes but it seems that a VMS with non-negligible 
metallicity, if left to evolve on its own, will lose n iost of its mass by s tellar winds 
and produce onlv a rather low-mass remnant (Belkus et al. 20071: Yungelson 



et al. I2OO7I ). The collisional growth can also be self-limited if the average time 
between two mergers becomes shorter than the thermal timescale of the VMS, 
as se en in Monte-Carlo simulations of verv rich cluster of single s tars fFreitag 
et al. l2006ar ). Then the VMS can not contract back to the MS and might instead 
become too diffuse to stop impactors. 

It is in any case unlikely that this scenario can lead to the formation of a 
continuous spectrum of massive stars. Simulations as well as simple mathemat- 
ical models suggest that one (possibly two) VMS detaches itself from the initial 
MF with very few, if any, stars populating the gap between mmax and rrir 



4. Second Scenario: Accretion-Driven Core Collapse 

The relaxation-driven scenario is based on a simplistic model of young clusters, 
neglecting the role of non-spherical substructure, initial mass segregation and 
interstellar gas. The last point is of particular importance because during the 
first Myr or so of their lives, most of the mass of clusters is in gas rather than 
stars. When the residual gas is expelled through ionisation by the new born 
massive stars or, for the most massive clusters, by the first SN explosion, the 
cluster reacts by expanding consider ably, probably to the point of complete 
dispersion for mos t low-mass c lusters ("Gever & Burkert"200l'; Kroupa, Aarseth, 
& Hurley 2001,; Goodwin k Bastia n 2006: Baumga rdt fc Kroupa 2007i ). If a 
cluster which remains bound is a factor Xji larger and a factor Xm less massive 
than in the embedded phase, with A/j ~ 10 and A^r ~ 0.3 being reasonable 
values, one can estimate that its relaxation time has increased by a factor Aj-ix ~ 

3 /2 

A^ with I' > 0.5 depending on what fraction of the stars are lost. Arix could 
be as large as ~ 20, indicating that the segregation-driven core collapse time 
could be much shorter in the embedded phase than estimated from observations 
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of young, gas- free clusters. Incidentally, the collision time increases by X^on ~ 
Xj^ X^,j with v' < 0.5, a factor which can reach hundreds. 

Such estimates suggest that the focus of the study of collisions in young 
clusters should be shifted to the embedded phase and to stars that are still 
in the process of formation, rather than MS objects. In this context, another 
process can lead to a strong increase of the stella ,r density, possibly up to a 
cohisional phase (|Bonnen. Bate, Zinnec kcr 19981: iBonnell Batd [20021 ) . The 



orbits of stars accreting from the interstellar gas shrink as a consequence of 
conservation of linear momentum. If the gas has zero net momentum, does 
not dominate the potential, and is accreted on a timescale long compared to 
orbital time, the size of the stellar system will contract like i?sys -^sysi where 
Mgys is the combined, increasing mass of the stars. The above assumptions 
should be valid for the core of a forming cluster or sub-cluster, rather than the 
whole cluster, A'gys < A^. The overall contraction of the system is eventually 
halted by small- stellar dynamical effects, in particular energy release by hard 
binaries. The maximum stellar density that can be reached is estimated to be 
Pmax = {n^)nina.x ~ {^'hot / {fn))^ P, where (m) is the mean stellar mass reached 
at that time in the contracting system, Mtot the, total mass of the parent proto- 
cluster and p its average mass density. This density appears to be high enough 
for collisions to affect a significant fraction of stars only in rather rich systems, 
Xsys ^ 100 (Clarke & Bonnell, in preparation). 

Clearly the role of collisions in massive star formation or, more generally, 
in the dynamise of young clusters can not be assessed independently of all other 
processes taking place in these complex systems. In this c ontribution, I have 
surveyed a few aspects that re quire consideration (see, e.g., iBallv &: Zinneckei 



I2OO5I : IZinnecker &: Yorkil2007l for other discussions of this subject). It seems 



that a complete picture will only be reached by means of numerical simulations 
combining stellar dynamics, the physics of the interstellar gas and the formation 
and early evolution of single and binary stars. 
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